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On a weak Gaufi law in general relativity and torsion 
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7'- Abstract 



We present an explicit example showing that the weak GauB law of general relativity (with 
cosmological constant) fails in Einstein-Cartan's theory. We take this as an indication 



^ ' that torsion might replace dark matter. 
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1 Introduction 

GauB' law in Maxwell's theory allows to link the total charge of an isolated charge distri- 
bution to the asymptotic strength of its electric field. In general relativity a weak form 
of Gaufi' law holds in the static and spherically symmetric case and remains true when a 
cosmological constant is added. We show that this weak GauB law fails when torsion is 
added a la Einstein-Cartan [H El [3]. 

2 The weak Gaufi law of general relativity 

To warm up, let us quickly review the derivation of the weak GauB law in general relativity 
with cosmological constant. 

In a first step we solve the Killing equation 

in the static, spherical case, i.e. for the four vector fields: ^ = d/dt generating time 
translation and 

(2) 

(3) 
(4) 

generating the rotations around the x-, y- and 2;-axis. In fact, since for the infinitesimal 
rotations, the commutator of the last two gives the first, we may discard the infinitesimal 
rotation around the x-axis. It is a typical accountant's work to list and solve the remaining 
3x10 linear, first order partial differential equations for the 10 unknowns g^v The solution 
is well known: 

dr^ = Bdt^ + 2D dtdr-A dr^ - C dO'^ - C sin^ 6 d(/?^ (5) 

with four functions A, 5, C and D of r. By a suitable coordinate transformation we may 
achieve D = and C = r^. Then A and B are positive. 
In a second step we solve the Einstein equation 

Ricci^^ - I scalar g^,^ - kg^u = SttG t^^, (6) 

in the vacuum, r^;^ = 0, and obtain the Kottler or Schwarzschild - de Sitter solution: 

B=-=l---\Kr\ (7) 

A r ^ ^ ' 

with an integration constant 5*. 
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In a third step we use the weak GauB law to determine the integration constant S in 
terms of the total mass of a spherically symmetric mass distribution T^t{f) with support 
inside a ball of radius R (not to be confused with the curvature scalar, that we denote by 
'scalar'). This law is simply the *i component of the Einstein equation: 

Ricci*i - i scalar g\ - Kg\ = 8nG t\ . (8) 

In the static, spherical case this equation reduces to: 



1 



1- 



d r 
dr A 



- A = SttG t' 



t ■ 



(9) 



Solving for d{r/A)/dr and integrating we have: 



= r-2G I r*t(f ) Anf^ df - |A r^ + K. 



(10) 



The integration constant K is seen to vanish by evaluating at r = and noting that ^4(0) 
is positive. Then the interior solution is 



A{r) 



1 - ( 2G I T\{r) 47rf2 drj /r - |A r^ 
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and from the continuity oi A at r = Rwe obtain the Schwarzschild radius S = 2GM with 
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(12) 



3 The weak Gaufi law in Einstein-Cartan's theory 

Our task is to redo the above three steps including torsion a la Einstein-Cartan. 



3.1 First step: invariant connection 

We have to solve the analogue of the Killing equation [4J for the now independent (metric) 
connection F 



d 



_ d^ 



ca ^ pA _ '^'^ pA _|_ '-"^ pA _|_ ^S -pA 



d^ 



T\.n + ,A4— = 0, 



(13) 



for the vector fields generating time translation and rotations. These 3 x 64 equations 
yield the following non-vanishing connection components: 



r^fec — x"'bc{r)-, 





aip 



a, b, c e {t,r}, 

n.fl = sineF"(r), 



smO cos 6, T'^dyy = T'^^^e = cos6'/sin^. 
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(14) 
(15) 
(16) 
(17) 



with 8 + 2 + 2 + 2 + 2 = 16 arbitrary functions X, E, F, C, Y oir. The metricity condition, 

d 



r. X 3iJ.v — ^^fixOfiu — ^"uxg^u — 0, (18) 

reduces the 16 functions to four arbitrary functions Ct, Cr, Dt, Dr and we have the fol- 
lowing non-vanishing connection components: 

r*,, = '^B'/B, r\a = D, A/B, T\, = \A/A, T\, = D^, (19) 

r'ee = rV^/ sin^ 6 = Q r^B, Vee = T\^/ sin^ 9 = -C, r^A, (20) 

r'ae = r\^ = Ca, T%r = T%r = 1/r, (21) 

T\^ = - sin ^ cos 9, V^e^ = ^'^^e = cos 0/sin 6. (22) 

3.2 Einstein's equation 

So far we have used a holonomic frame dx^. Now it will be convenient to work in an 
orthornormal frame e" =: e'^^dx'^. We will use the notations of reference [5]. The metric 
tensor reads gf,u{x) = e^^i^x) e^(a;) rjab- 

The connection with respect to a holonomic frame is written as a 5'^(4)-valued 1-form 
r°^ =: r^^^dx'^. The link between the components of the same connection with respect 
to the holonomic frame F and with respect to the orthonormal frame u is given by the 
GL(4) gauge transformation with e(x) = e"^(x) G GL{4:); 

uj = eTe-^ + ede-\ (23) 

or with indices: 

dxf^ 

For u the metricity condition is algebraic and means that its values uj°'b are in the Lie 
algebra of the Lorentz group: Uab = —^ba- 
in, the orthonormal frame with e = diag(v5, v^, ?", rsin^), the non- vanishing com- 
ponents u:°'biJL of the invariant connection flT9l - l22|) are: 

u\t = Dt^A/B, uj\.r = Dr^/A/B, u\e = uj\^/sine = Ctr/^, (25) 

w^w = ~ cos 6, uj^'ee = ^'' wl sin d = -G,- r/y/A. (26) 

The curvature as so(l, 3)-valued 2-form, 

R:=du + ^[uj,u], (27) 

has the following non-vanishing components, W^b ='■ ^R'^bfj.udx'^dx'^: 

R\tr = -{Dt^/A/B)', R'ete = -C.A t/VB, R\.,e = (Q r/v^)' - C,.D, rjyfB, 

R iptip = sin 6 R 0te, R ipnp = sin 6 R ere, 
Wete = CtDt r^A/B, R'ore = -(C. r/VI)' + CtDr r^/B, 
R^ifitifi = ^^^0 R^dtd, R^^rifi = s^^S R^ere, (28) 

R%e^ = sin^ (1 + C^r^B - C^.r^A). 



oo\^ = e\ T-p, e-'^ + e\ -^^ e'^^. (24) 



Next we compute the Ricci tensor, 

Ricci^ := ri^^'R^'^.e-^^ce-^^, 

whose non-vanishing components are: 

Ricci*t = (A \/A/B)'/^/AB + 2 CrDt/B, 

Ricci", = {Dt^/A/By/VAB + 2{rCr/VAy/{rVA)-2CtDr/B, 

RicciV = -2{Ct/^/By/^/A-2Ct/{r^/AB) + 2CrDr/^/AB, 

Ricci'^t = -2CtDty^A/B/B, 



Ricci^ 



Ricci"^,, 



(29) 



(30) 
(31) 
(32) 
(33) 
(34) 



(a/v^)7VI - ly + Cr/{rA) - Cl/B + Cl/A - CtDr/B + CrDt/B. 

The curvature scalar is 

scalar = Ricci*( + Ricci", + 2 Ricci^e 

= 2 (A ^/AjByi^fAB + 4 (a/ VI)7v^ - 2/r2 + 4 ^/(rA) 

-2 C^/fi + 2 C^/A - 4 a A/5 + 4 a A/5. (35) 

In an orthonormal frame the Einstein equation reads 



i?"b - \ scalar 5% - A 5% = SttG r^ 



(36) 



In our orthonormal frame e = diag(v5, V^ 
is symmetric and 



, r, rsin^), the energy momentum tensor 



r„ 



T b 



I p{r) q{r) \ 

— g(r) — Pr(^) 

-pa(r) 

V -Pa(r)/ 

The tt, rr and ^^ components of the Einstein equation are: 

- 2 (a/VI)7 VI + l/r^ - 2 a/(rA) 

+C,V5 - C,VA + 2 Ct A/5 - A = SttGp, 
l/r2 + QVfi-C,VA-2aA/5-A = -SttGp,, 
-(A ^'A/By/^AB - {Crl^fAyi^fA - Cr/{rA) 

+CtDr/B -CrDt/B -K = -SyrGpa- 

The two off-diagonal components read 



(37) 
(38) 



(39) 
(40) 

(41) 



- 2 {Ct/VBy/VA - 2 Gt/(r\/AB) + 2 G^ A/v^AB = 2 QA ^AjBjB = SnC q. (42) 



3.3 Cartan's equation 

The torsion as M^-valued 2-forni, 

T := Be = de + ue, (43) 

has the following non-vanishing components, T" =: Ta'bcV'^ "e^e'^, 

Tttr = (A - IB'/A) VA/B, Trtr = Dr/Vs, (44) 

Tete = T^t^ = Ct/^fB, Tore = T^r^ = (C, - l/r)/v^. (45) 

The Cartan equation, 

T'^e'^eabcd = -SnG Sab, (46) 

determines the torsion in terms of its source, the half- integer spin current. This is the 
Lorentz-valued 3-form Sab i-e. the variation of the matter Lagrangian with respect to the 
spin connection Uab- To simplify the Cartan equation, let us decompose the torsion tensor 
into its three irreducible parts: 

Tabc = Aabc + VabVc - llaM + Mabc, (47) 

with the completely antisymmetric part Aabc '■= |(^afec + Tcab + Tbca), the vector part 
Vc ■■= iTabcV"^, and the mixed part Mabc characterized by Mabc = -Macb, MabcV"^ = 0, 
and Mabc + Mcab + Mbca = 0. Likewise, we decompose the spin tensor Sabc defined by 

Sabc = Clabc + VcaSb " VcbSa + ^riabc, (48) 

with the completely antisymmetric part aabc '■= ^{sabc + Scab + Sbca), the vector part 
Sb '■= ^SabcV"''^, and the mixed part rriabc characterized by rriabc = —f^bac, i^abcV"''^ = 0, and 

mabc + mcab + mbca = 0. 

Then the Cartan equation reads: 

Aabc = -8nG ttabc, Va = ISnGSa, Mcab = -^T^G 171 abc- (49) 

In the static, spherical case, we have Aabc = 0, 

Vt = l{2Ct + Dr)/VB, Vr = \{-\B'/B-2/r + 2Gr + DtA/B)/^, (50) 

and 

Mrtr = i(-a + Dr)l^, Mt^t = U\B' / B - l/v + C, " DtA/B)/VA, (51) 
Mete = M^t^ = -\Mrtr, Mere = M^,^ = ^Mtrt- (52) 

3.4 Second step: vacuum solution 

In vacuum. Tab = and Sabc = 0, we retrieve the Kottler solution, equation ((Tj). Indeed 
by Cartan's equation, vanishing spin current implies vanishing torsion: Ct = D,. = 0, 
Gr = 1/t and Dt = jB'/A and then the invariant, metric connection (IT9]- [22l) reduces to 
the (symmetric) Christoffel symbols. 
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4 Third step: a Schwarzschild star with torsion 

Note that metric and connection in the static, spherical case are automatically invariant 
under space inversion. (This is not true in the homogeneous, isotropic case [HE].) However 
we do not have invariance under time reversal and to be precise we should say 'stationary' 
rather than 'static'. 

To construct a counter example to the weak GauB law, we do suppose invariance 
under time reversal. This is certainly not justified for our sun, but not unreasonable 
for a Schwarzschild star with constant mass density dp/dr = inside the radius R. 
The following functions are odd under time reversal and must vanish: all connection, 
curvature, Ricci, energy-momentum, torsion and spin tensor components with an odd 
number of indices equal to t. Consequently Ct, Dr, q, St and m^r are zero and we 
remain with four unknown functions of r in the fields: B, A, C^ and Dt. In the sources 
we still have five arbitrary functions of r: the mass density p the radial and azimuthal 
pressure pr, Pa and the spin densities Sr and mtrt- They define the right-hand sides of the 
five remaining field equations, three Einstein and two Cartan equations. 

To continue, we set Dt = \B' /A and simplify notations C := C^, s := Sr ■ Then the 
two Cartan equations reduce to: 

2(C- l/r)/yi=87rGs, mtrt = \s. (53) 

Now we may introduce a Schwarzschild star by assuming that the mass density p and the 
spin density s are constant, i.e. r- independent, with an equation of state: s = wp. 

Upon eliminating C via the Cartan equation, the tt, rr and 66 components fl39l - HTl) 
of Einstein's equation reduce to: 



A'/irA^) + l/r^ -l/{r^A) -lQ7rGwp/{rVA) - {AnGwp)'^ - A = 8nG p, (54) 

l/{rA) + 47TGwp/y/A] B'/B 

-1/r^ + l/{r^A) +8nGwp/{r^/A) + {AttGwpY + A = SnGpr, (55) 

'^B"/{AB) - \iA'/A + B'/B) B'/{AB) 

-UA'/A-B'/B)/{rA) + aB'/B + l/r)AnGwp/y/A + A = SnGp^. (56) 



As with zero torsion, the tt component decouples from the other two equations and can 
be integrated separately. To redo the third step of section 2 for the tt component with 
torsion, we now need two definitions of mass: an interior mass, 

rR 

Mi := / p 47rf 2 df = InR^ p, p = T^'=\=t = r''=\=u (57) 

and an exterior mass Mg defined by the strength of the gravitational field outside, r > R, 



A{r) 






(5^ 



In contrast to the torsionless case, they do not coincide: 



Mp = M 



Qw 



R 



fdr 

7w) 



+ 



3^2 CM 



2i?3 



(59) 



Note that for sufficiently large \w\ the exterior mass exceeds the interior one, even for 
negative w. Note also that this mass relation depends on the interior solution A{r) of 
the tt component of the Einstein equation (15^ . This solution is not obvious (to us) and 
we will solve equation (1541) numerically. As a test the numerical solution will reproduce 
for w; = the Schwarzschild star with cosmological constant [TJ [U [9]. This solution has 
Pr = Pa. ='■ P- Its functions. 



1 



A=—, B = {aK + (5W)\ p = p 



K 



are continuous at the boundary r = R. The auxiliary quantities used are: 
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^TiGp + A), a:=i87rGp/7, /3 := (-iSvrGp + iA)/7 = 1 



W{r) := ^l--fr\ K := W{R). 



a, 



(60) 



(61) 
(62) 



5 Numerical solution 

Equation flS^ has an integrable singularity at r 
the dependent variable, a{r) := r/A{r), yielding: 



which can be avoided by redefining 



IQTiGwpy/ra — {AirGwprY — Ar^ = SnGpr'^. 



(63) 



We solve this equation by a Runge-Kutta algorithm with initial condition a(0) = for 
< r < R. We check that a{r) remains positive for < r < R, and that linir^o A{r) = 1 
which also ensures that that r/ ^^/A(r) = ^r a{r) is integrable for < r < i?. Then we 
get the masses from 



Mi = ^nR^p, Me = {R- \AR^ - a(i?))/(2 G) . 



(64) 



For w = we reproduce the analytic solution A = [1 — |(87rGp + A)r^]~^ and have 
Mi = Me in accordance with the weak GauB law. To obtain a ratio of Me/Mi = 5 for 
the sun, M, = Mq, R = 7 ■ 10^ m, we must choose w = 3.1 s. We get the same ratio 
for a cluster, Mj = W^^ Mq, R = 3 ■ 10^^ m, with w = 2 ■ 10-*^^ s. In these two cases 
the positive definite contribution ^w"^ GMi/R^ to the mass ratio amounts to 6 • 10^^ and 
0.6 respectively. We have used the experimentally favoured value of A = 1.5 ■ 10~^^ m~^. 
Setting the cosmological constant to zero however does not change the values of Mg/Mj 
by more than 10~^. 



6 Conclusion 

A torsion induced failure of the weak GauB law might be welcome with respect to the 
dark matter problems. Indeed we have already seen [1] that the Hubble diagram of super 
novae can be fitted by the Einstein-Cartan theory with w = 10^^ s and no dark matter. 
This w-value is not far from the one found here for a spherical cluster. However they are 
far, far away from the naive microscopic value: 

w = — - ~ 10~2^ s. (65) 

It would nevertheless be interesting to compute the rotation curve of a realistic galaxy 
and lensing in the Einstein-Cartan theory. Both are formidable theoretical challenges. 
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